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Abstract 

We discuss a necessary nonvalence contribution in time- 
like exclusive processes. Following a Schwinger-Dyson 
type of approach, we relate the nonvalence contribution to 
an ordinary light-front wave function that has been exten- 
sively tested in the spacelike exclusive processes. A com- 
plicate multi-body energy denominator is exactly cancelled 
in summing the light-front time-ordered amplitudes. Ap- 
plying our method to K ^3 and D° — > K~~t + vi where a 
rather substantial nonvalence contribution is expected, we 
find not only an improvement in comparing with the exper- 
imental data but also a covariance(i.e. frame-indepndence) 
of existing light-front constituent quark model. 

1 INTRODUCTION 

As discussed in this PEP-N meeting, the facilities that co- 
piously produce the lower-lying mesons such as K and D 
can provide a lot of rich physics as exciting as the new 
and upgraded _B-meson factories promise. To fulfill such 
excitement, however, intensive theoretical studies should 
be accompanied in the analyses of exclusive meson decays 
and form factors. Thus, more and more scrutinized model 
analyses are called for. 

Perhaps, one of the most popular formulations for the 
analysis of exclusive processes may be provided in the 
framework of light-front (LF) quantization [ffh. In partic- 
ular, the Drell-Yan-West (q + = q° + q 3 = 0) frame has 
been extensively used in the calculation of various elec- 
troweak form factors and decay processes [^, [| |], 0], As 
an example, only the parton-number-conserving (valence) 
Fock state contribution is needed in q + =0 frame when the 
"good" component of the current, J + or Jj_ = (J x , J y ), 
is used for the spacelike electromagnetic form factor cal- 
culation of pseudoscalar mesons. The LF approach may 
also provide a bridge between the two fundamentally dif- 
ferent pictures of hadronic matter, i.e. the constituent quark 
model (CQM) (or the quark parton model) closely related 
to the experimental observations and the quantum chromo- 
dynamics (QCD) based on a covariant non-abelian quan- 
tum field theory. The crux of possible connection between 
the two pictures is the rational energy-momentum disper- 
sion relation that leads to a relatively simple vacuum struc- 
ture. There is no spontaneous creation of massive fermions 
in the LF quantized vacuum. Thus, one can immediately 
obtain a constituent-type picture, in which all partons in 



a hadronic state are connected directly to the hadron in- 
stead of being simply disconnected excitations (or vacuum 
fluctuations) in a complicated medium. A possible realiza- 
tion of chiral symmetry breaking in the LF vacuum has also 
been discussed in the literature [Q]. 

On the other hand, the analysis of timelike exclusive pro- 
cesses (or timelike q 2 > region of bound-state form fac- 
tors) remained as a rather significant challenge in the LF 
approach. In principle, the q + ^ frame can be used 
to compute the timelike processes but then it is inevitable 
to encounter the particle-number-nonconserving Fock state 
(or nonvalence) contribution. The main source of difficulty 
in CQM phenomenology is the lack of information on the 
non-wave-function vertex(black blob in Fig. |j](a)) in the 
nonvalence diagram arising from the quark-antiquark pair 
creation/annihilation. The non-wave-function vertex(black 
blob) was recently also called the embedded state [^]. This 
should contrast with the white blob representing the usual 
LF valence wave function. 
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Figure 1 : Effective treatment of the LF nonvalence ampli- 
tude. 

In principle, there is a systematic program laid out by 
Brodsky and Hwang []8|] to include the particle-number- 
nonconserving amplitude to take into account the nonva- 
lence contributions. However, such program requires to 
find all the higher Fock-state wave functions while there 
has been relatively little progress in computing the basic 
wave functions of hadrons from first principles. Recently, 
a method of analytic continuation from the spacelike region 
has also been suggested to generate necessary informations 
in the timelike region without encountering a direct cal- 
culation of the nonvalence contribution [|]]. Even though 
some explicit example has been presented for manifestly 



covariant theoretical models, this method has not yet been 
implemented to more realistic phenomenological models. 

In this talk, we thus present an alternative way of han- 
dling the nonvalence contribution. Our aim of new treat- 
ment [ pi] ] is to make the program more suitable for the 
CQM phenomenology specific to the low momentum trans- 
fer processes. Incidentally, the light-to-light (K^) and 
heavy-to-light (D° —> K~£ + v$) decays involving rather 
low momentum transfers bear a substantial contribution 
from the nonvalence part and their experimental data are 
better known than other semileptonic processes with large 
momentum transfers. Including the nonvalence contribu- 
tion, our results on Ka and D° — > K~l + vi not only show 
a definite improvement in comparison with experimental 
data but also exhibit a covariance (i.e frame-independence) 
of our approach. 

This talk is organized as follows. In Section 2, we 
present the non-wave-function vertex in the nonvalence di- 
agram in terms of light-front vertex functions, utilizing the 
covariant Bethe-Salpeter(BS) model of (3+l)-dimensional 
fermion field theory. The nonvalence part of the weak form 
factors for 0~ — > 0~ semileptonic decays is expressed 
in terms of light-front vertex functions of a hadron and a 
gauge boson. The link operator connecting (n — l)-body to 
(n+l)-body in a Fock state representation is obtained by an 
analytic continuation of the usual BS amplitude. We also 
show that the complicated (n + 2)-body energy denomina- 
tors are exactly cancelled in summing the light-front time- 
ordered diagrams. In Section 3, we show our numerical 
results for K(z and D° — » K~l + v? decays. Conclusions 
follow in Section 4. 

2 NEW EFFECTIVE TREATMENT 

2.1 0~ — ► 0~ Semileptonic Decays 

The semileptonic decay of Qiq bound state with four- 
momentum and mass Mi into another Q2 q bound state 
with P-2 and M2 is governed by the weak current, viz., 

J M (0) - (P 2 |Q 2 7 M Qil^i} 

= U(q 2 )(P 1+ P 2 r + .Uq 2 )q^, (1) 

where q^ = (Pi — P2) M is the four- momentum transfer to 
the lepton pair (Iv) and m\ < q 2 < (Mi - M 2 ) 2 . The 
covariant three-point Bethe-Salpeter (BS) amplitude of the 
total current J M (0) in Eq. (jj]) may be given by 



J"(0) = tN c 



d*k 



(2TtY (jpl - ml + ie)(p 
1 



(2) 



where N c is the color factor, -ffj™ is the covariant inti- 
tial[final] state meson-quark vertex function that satisfies 
the BS equation, and = Tr[7 5 (^ + 7711)7^(^3 + 



As shown in the literature [Qj, the LF energy integra- 
tion reveals an explicit correspondence between the sum of 
LF time-ordered amplitudes and the original covariant am- 
plitude. For instance, performing the k~ pole integration, 
we obtain the LF currents, Jy and J^ v corresponding to 
the usual LF valence diagram and the nonvalence diagram 
shown in Fig. |l|(a), respectively. Since HZf™ satisfies the 
BS equation, we iterate H™** once and perform its LF en- 
ergy integration to find the corresponding LF time-ordered 
diagrams Figs. [j](b) and |l](c) after the iteration. The sim- 
ilar idea of iteration in a Schwinger-Dyson (SD) type of 
approach was presented in Ref. jlOt ] to pin down the LF 
bound-state equation starting from the covariant BS equa- 
tion. 

Comparing the LF time-ordered expansions before and 
after the iteration, we realize that the following link be- 
tween the non-wave-function vertex (black blob) and the 
ordinary LF wave function (white blob) as shown in Fig. [l] 
naturally arises, i.e., 

(M 2 - MlW ( Xl ^) 

[dy][d 2 \±]K(xi,]g.±i\ yj,\±j)^(y j ,\±j), (3) 



where M is the mass of outgoing meson and Ai 2 , = (m\ + 
kj_i)/ x i — ( m i + ki 2 )/(~ x 2) with xi = 1 — X2 > 1 due 
to the kinematics of the non-wave-function vertex. 
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given by p x = P x — k, p 2 = P2 — k, and p q = k. 



Figure 2: Non-wave-function vertex(black blob) linked to 
an ordinary LF wave function(white blob). 

We note that Eq. @) essentially takes the same form as 
the LF bound-state equation (similar to the LF projection 
of BS equation) except the difference in kinematics(e.g. 
— X2 > for the non-wave-function vertex). Incidentally, 
Einhorn jl2| ] also discussed the extension of the LF BS am- 
plitude in 1+1 QCD to a non-wave-function vertex similar 
to what we obtained in this work. 

In the above procedure, we also find that the four-body 
energy denominator D4 is exactly cancelled in the sum 
of LF time-ordered amplitudes as shown in Figs. |l](b) 
and 0(c), i.e., l/D 4 Z>f + l/D A D% = l/D%D%. We thus 
obtain the amplitude identical to the nonvalence contri- 
bution in terms of ordinary LF wave functions of gauge 
boson(W0 and hadron (white blob) as drawn in Fig. [jjd). 
This method, however, requires to have some relevant op- 
erator depicted as the black square(/C) in Fig. ||(See also 
Fig. |l](d)), that is in general dependent on the involved 
momenta connecting one-body to three-body sector. We 
now present some details of kinematics in the semileptonic 



decay processes to discuss a reasoning of how we han- 
dle the nonvalence contribution involving the momentum- 
dependent JC for relatively small momentum transfer pro- 
cesses such as 7r e 3, Kgz and D — > K Iv. 

2.2 Kinematics and Model Description 

Our calculation is performed in purely longitudinal mo- 
mentum frame jTT| , pj[ ] where q + > and Pu_ = P 2 ^ = 



so that the momentum transfer square q 



q + q~ > 



is time-like. One can then easily obtain q 2 in terms of 
the momentum fraction a = P 2 + 1 ' P\ = 1 — Q + / Pi~ as 
g 2 = (1 — a){Mi — Mf/a). Accordingly, the two solu- 
tions for a are given by 



Mi 



Mf + M| - 
2MiM 2 



'Mf + Mf -g 2 
2MiM 2 



(4) 

The +(— ) sign in Eq. (Q) corresponds to the daughter me- 
son recoiling in the positive(negative) z-direction relative 
to the parent meson. At zero recoil(<7 2 = g 2 lax ) and maxi- 
mum recoil^ 2 = 0), a± are given by 



a+^max) = a-(<2max) = 



a+(0) = l, a_(0) = 



M2 
Mi 

M2 
Mi 



(5) 



In order to obtain the form factors f±(q 2 ) which are in- 
dependent of a±, defining </ + (0)| Q=ct± = 2P 1 + H + (a±) 
from Eq. (jj]), we obtain 



M 2 



± 



(lTa+)ff + («-) 



(6) 

The form factors f+(q 2 ) and f-(q 2 ) are related to the 
scalar form factor fo(q 2 ) in the following way: 



/ + (0) = / (0), jo(<z 2 H/+(g 2 ) + 



Mf - M 2 



(7) 

The differential decay rate for 0~ — > _ semileptonic de- 
cay is given by [0] 



l \2 



X l [ [K f (q 2 )] 2 (l + ^)\f + (q 2 f 
M 9 2 9 3 m 2 . 

+m ^ 1 -m|)W /o(9)I 



(8) 



where is the Fermi constant, V qi g 2 is the element 
of the Cabbibo-Kobayashi-Maskawa(CKM) mixing matrix 
and the factor Kf(q 2 ) is given by 



K f (q 2 



1 



2Mj 



nl/2 



(M 2 + M 2 - q 2 ) 2 - 4M 2 M 2 



(9) 



With the iteration procedure Eq. (|3|) in this q + > 
frame, the results for the "+"-component of the current J p 
in Eq. (^) are given by 

^(z,k ± )S+%(z',k ± ) 



and 



16tt 3 Jo 



da; / d 2 k 



x(l - x)(l - x') 



(10) 
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x / 7~~ r / d 2 lx^(a;,k ± ;y,U)%(y,l ± )j(ll) 
J 2/(1 ~ 2/) 7 

where Eq. ([|) has been used for the nonvalence wave func- 
tion at the black blob as shown in Fig. f|(see also Fig. 0(d)). 

The ordinary LF vertex functions(white blob in Fig. [I]) 
in Eqs. ( JTo| ) and ( p"T| ) are given by 
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1 — x 
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k 2 m| + k 2 

-i ■- g ; \ (12) 



where x = k + / 'P± ', x' = x/a. The ty g in Eq. ( |l l| ) cor- 
responds to the light-front energy denominator(i.e. D s in 
Fig. 0(d)) and its explicit form is given by 

* s (z,k ± ) = — t— L- — . (13) 



l-a 



1-x 



We call ^ g the light-front vertex function of a gauge bo- 
son []. 

In Eqs. (|l0|) and (|l]), the trace terms Sy{p g = 
k m ) = (4P 1 +/a; / ){ki + [a^m + (1 - a;)m g -][i'm 2 + (1 - 
x'jmg]} and S^ v jpJ" = pi~on) = s v(p7 = P»J + 
^lojjP^nfe; — -P9011) correspond to the product of initial 
and final LF spin-orbit wave functions that are uniquely de- 
termined by a generalized off-energy-shell Melosh trans- 
formation. Here, the subscript (on) means on-mass-shell 
and the instantaneous part of nonvalence diagram corre- 
sponds to 4p+ onP + on ( P q -pq on ) in While the LF ver- 
tex function ft,fr|j formally stems from H^Sj, practical in- 
formations on the radial wave function ^jr/i (x, kj_) (con- 
sequently hfn-i) can be obtained by LF CQM. The details 
of our variational procedure to determine both mass spectra 
and wave functions of pseudoscalar mesons were recently 
documented in Refs. [g, 0] along with an extensive test of 
the model in the spacelike exclusive processes. The same 
model is used in this work, i.e., comparing the LF vertex 
functions ^ in Eq. (|l2|) with our light-front wave function 
given by Ref. [|[ Q], we identify 



N' 



0:r 



M n 



(14) 

1 While one can in principle also consider the BS amplitude for \& 9 , 
we note that such extension does not alter our results within our approx- 
imation in this work because both hadron and gauge boson should share 
the same kernel. 



rriq) 2 and the Jacobian of 



where M = y-Mg — (m q 
the variable tranformation k = (k z , kj_) — > (a;, kj_) is ob- 
tained as dk z /dx = Mq/\Ax{\ — x)} and the radial wave 
function is given by 



3 NUMERICAL RESULTS 
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exp(-k 2 /2/3 2 ), (15) 



which is normalized as J d 3 fc[0(k 2 )| 2 = 1. Substituting 
Eqs. (0) and ([[J) into Eqs. @ and (0), one can obtain 
the valence and nonvalence contributions to the weak form 
factors for — > CP semileptonic decays in light-front 
quark model. 

While the relevant operator K, is in general dependent 
on all internal momenta (x, k^, y, a sort of average on 
K. over y and 1^ in Eq. (11) which we define as Gp 1 p 2 = 



f[dy] [d 2 l±]K.(x, kj_; y, l±)"i>f(y, is dependent only on 
x and k^ . Now, the range of the momentum fraction x de- 
pends on the external momenta for the embedded states. 
As shown in Eq. ([[l]), the lower bound of x for the ker- 
nel in the nonvalence contribution is given by a which has 
the value a = M^jMx at the maximum q 2 . As the mass 
difference between the primary and secondary mesons gets 
smaller, not only the range of q 2 is reduced but also a gets 
closer to 1 . Perhaps, the best experimental process for such 
limit may be the pion beta decay tt^ 1 — ► ir°e ± v e , where our 
numerical prediction /_(0)//+(0) = -3.2 x 1(T 3 follow- 
ing the treatment presented in this work is in an excellent 
agreement with —3.5 x 10 -3 obtained by the method pro- 
posed by Jaus [|l4j] including the zero-modes JT6[ , [| O], 
In Ademollo-Gatto's SU(3) limit [pj, the q 2 range of the 
nonvalence contribution shrinks to zero and a becomes pre- 
cisely 1. However, even if a is not so close to 1, the initial 
wavefunction &i(x, kj_) plays the role of a weighting fac- 
tor in the nonvalence contribution and enfeeble the contri- 
bution from the region of x near 1 . Thus, for the processes 
that we discuss in this talk, the effective x region for the 
nonvalence contribution is quite narrow. Similarly, the re- 
gion of the transverse momentum k^ is also limited only 
up to the scale of hadron size due to the same weighting 
factor kj_). Here, we thus approximate Gp 1 p 2 as a 
constant and examine the validity of this approximation by 
checking the frame independence of our numerical results. 

For the check of frame-independence, we also compute 
the "+" component of the current in the Drell-Yan- 
West (q + = 0) frame where only valence contribution 
exists. Since the form factor f+(q 2 ) obtained from 
in q + — frame is immune to the zero-mode contribu- 
tion || U m H], the comparison of f+(q 2 ) in the two 
completely different frames (i.e. q + — and q + ^ 0) 
would reveal the validity of existing model with respect to 
a covariance. The comparison of f_(q 2 ), however, cannot 
give a meaningful test of covariance because of the zero- 
mode complication as noted in Ref. [Q. Indeed, the dif- 
ference between the two (q + = and q + ^ 0) results of 
f-(q 2 ) amounts to the zero-mode contribution. 
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Figure 3: The weak form factors for K® 3 compared with 
the experimental data fllTy, 

In our numerical calculation for the processes of 
and D° — » K~l + v decays, we use the linear potential 
parameters presented in Ref. fl]. In Fig. 0, we show the 
weak form factors f+{q 2 ) and fo(q 2 ) for K® 3 decays. The 
thick solid lines are our analytic solutions obtained from 
the q + = frame; note here again that the lower thick 
solid line (/o) in Fig. ||is only the partial result without in- 
cluding the zero-mode contribution while the upper thick 
solid line (/ + immune to the zero-mode) is the full re- 
sult. The thin solid lines are the full results of our ef- 
fective calculations with a constant (0^=3.95) fixed by 
the normalization of / + at q 2 = limit. For comparison, 
we also show only the valence contributions(dotted lines) 
in q + 7^ frame. As expected, a clearly distinguishable 
nonvalence contribution is found. Following the popular 
linear parametrization [jjJJ], we plot the results of our ef- 
fective solutions(thin solid lines) using fi{q 2 ) = fi{q 2 = 
m 2 )(l + Xiq 2 /M 2 + )(i — +,0). In comparison with the 
data, the same normalization as the data /+(0) = 1 Jl7| ] 
was used in Fig. ||. Our effective solution(upper thin solid 
line) is not only in a good agreement with the data Jl7| ] 
but also almost identical to that in q + = frame(upper 
thick solid line) indicating the frame-independence of our 
model. Note also that the difference in fo(q 2 ) between 
q + ^ 0(lower thin solid line) and q + — 0(lower thick solid 
line) frames amounts to the zero-mode contribution. 

In comparison with experimental data, we summarized 
our results of several experimental observables in Table 1 ; 
i.e. the actual value of /+(0), the slopes A + [Aq] of f+(q 2 ) 



Table 1: Model predictions for the parameters of K® 3 decays. The decay width is in units of 10 6 s l . The used CKM 
matrix is \V„J = 0.2196 ± 0.0023 [ll 





Effective 


+ n 
= 


Experiment 


/+(0) 


0.962 [0.962] 


0.962 [0.962] 




A+ 


0.026 [0.083] 


0.026 [0.026] 


0.0288 ± 0.0015 [K^] 


A 


0.025 [-0.017] 


0.001 [-0.009] 


0.025 ± 0.006 [Jf" 3 ] 


U 


-0.013 [-1.10] 


-0.29[-0.41] 


-0.11±0.09[A^ 3 ] 




7.3 ±0.15 


7.3 ±0.15 


7.5±0.08 


r(^ 3 ) 


4.92 ±0.10 


4.66 ±0.10 


5.25±0.07 



[fo(q 2 )] at q 2 = 0, U=f-(°)/f+(°)> and the decay rates 
T(K%) and T(K° 3 ). In the second column of Table 1, 
our full results including nonvalence contributions are pre- 
sented along with the valence contributions in the square 
brackets. In the third column of Table 1, the results in 
q + = frame are presented with[without] the instanta- 
neous part. As one can see in Table 1, adding the nonva- 
lence contributions clearly improves the results of Ao, i.e. 
our full result of Ao= 0.025 is in an excellent agreement 
with the data, Ag Xp = 0.025±0.006. Since the lepton mass 
is small except in the case of the r lepton, one may safely 
neglect the lepton mass in the decay rate calculation of the 
heavy-to-heavy and heavy-to-light transitions. However, as 
one can see from the improved result for decay rate, 
the reliable calculation of fo{q 2 ) is required especially for 
since the muon(/i) mass is not negligible, even though 
the contribution of fo(q 2 ) is negligible for K e3 case. 

In Figs. 4(a,b), we show the weak form factors for 
D° — > K~ ' l + v decays and compare with the experimental 
data [|l8|] (full dot) with an error bar at Q 2 — as well as the 
lattice QCD results JT^] (circle and square) and [p0|](cross). 
All the line assignments are same as in Fig. 3. In Fig. 4(a), 
the thin solid line of our full result in q + ^ is not visi- 
ble because it exactly coincides with the thick solid line of 
the result in q + = confirming the frame-independence 
of our calculations. Our value of / + (0)=0.736 is also 
within the error bar of the data [Jl|], /+ xp (0)=0.7±0.1. 
In Fig. 4(b), the difference between the thin and thick 
solid lines is the measure of the zero-mode contribution 
to fo(q 2 ) in q + = frame. The form factors obtained 
from our effective calculations(Gi)^=3.5) are also plotted 
with the usual parametrization of pole dominance model, 
i.e. f +(Q) (q 2 ) = f +(Q) (0)/(l-q 2 /Ml (0+) ). Our pole 
masses turn out to be A/ 1 -=2.16 GeV and AT +=2.79 GeV, 
respectively, and we note that M 1 -=2.16 GeV is in a good 
agreement with the mass of £>*,i.e. 2.1 GeV. Using CKM 
matrix element | V cs \ = 1 .04 ± 0. 16 [180, our branching ra- 
tios Br(£>° 3 ) = 3.73 ± 1.24 and Br(L>° 3 ) = 3.60 ± 1.19 
are also comparable with the experimental data3.64±0.18 
and 3.22 ± 0.17 respectively. 

In Fig. U we show the differential decay rates for D° — > 
K ~e + v e and D° — ► K~ y^v^ transitions obtained from 
our effective solutions. As in the case of Kg 3 decays, we 
were able to evaluate the fo(q 2 ) contribution to the total 
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Figure 4: The weak form factors for D - 
The same line codes are used as in Fig. 3. 
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manner although its contribution is more suppressed than 
the K^z case. 

4 CONCLUSION 

In summary, we presented an effective treatment of the LF 
nonvalence contributions crucial in the timelike exclusive 
processes. Using a SD-type approach and summing the 
LF time-ordered amplitudes, we obtained the nonvalence 
contributions in terms of ordinary LF wavefunctions of 
gauge boson and hadron that have been extensively tested 
in the spacelike exclusive processes. Including the non- 
valence contribution, our results show a definite improve- 
ment in comparison with experimental data on K#3 and 
D° — > K~l + Vt decays. Our result on 7r e 3 is also consis- 
tent with the result obtained by other method. Furthermore, 
the frame-independence of our results indicate that a con- 
stant Gp 1 p 2 is an approximation appropriate to the small 
momentum transfer processes. A similar conclusion was 
drawn in a recent application of our method to the skewed 
quark distributions of the pion at small momentum transfer 
region pT|]. Applications to the heavy-to-light decay pro- 
cesses involving large momentum transfers would require 
an improvement on this approximation perhaps guided by 
the perturbative QCD approach. Consideration along this 
line is underway. 
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